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I . INTRODUCTION 


The physical and mathematical bases of the kinetic theory approach to 
the chemically reacting, incompressible, turbulent flows have been described 
in several of the previous papers.^ ^ Recently, it has been shown^^ that 
certain exact solutions of the governing kinetic equations can be obtained 
through construction of the appropriate Green's functions. However, the com- 
plexity of the kinetic equations is such that it is generally impractical to 
attempt to obtain the exact solutions except for the rather simple flow con- 
figurations such as the plane, shear mixing layer analyzed in Reference 6. 

A study is being conducted to develop ar. approximate formalism which 
embodies the salient contributions of the kinetic theory but which employs 
a much simplified set of equations such that the chemical reactions taking 
place in the realistic flow fields can be analyzed. 

This need to generate an approximate method of application of the kinetic 
theory to engineering problems becomes accentuated when one considers the 
compressible kinetic theory which is being constructed. Obviously, the 
compressibility increases the complexity of the governing equations, and, 
exact solution of the equations for the realistic flow configurations 
becomes extremely difficult. 

The essence of development of the approximate method of solution is as 
follows. 

We first develop a simplified moment relationship, from the two non- 
equilibrium-degree kinetic equation, ^ which will relate the effective 
length scale, to the mean turbulence energy, and perhaps a few other 
mean flow properties, and to the two given length scales, A^ and A^. This 
length-scale equation is not needed for the simple flow fields such as 
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the plane shear layer. 

Next, we make use of the existing solution of the momentum equation 
for the particular flow problem. This is because the greatest difficulty 
in the application of the kinetic theory arises from the solution of the 
governing equation for the turbulent momentum field when the flow field is 
complicated. At the same time, it is known that the mean velocity, mean 
shear stress, etc., are usually quite well described by the existing 
theories such as the Prandtl’s mixing-length theory. We assume an appro- 

*> -h 

priate form of the distribution function, f(t,x,u)du, and from the existing 
solutions of the mean values, we determine the distribution function of 
the fluid elements. 

Finally, the kinetic equations governing the chemical species and 

energy are solved for their distribution functions by the use of the 

f(t,x,u)du, since it is the mixing and chemical reactions which the exist- 

( 1 - 3 ) 

ing theories are unable to describe correctly. The effective length- 

scale equation is also employed in the solution of the kinetic equations 
when it is warranted. 

The above process is simpler when the incompressible assumption is 
made, for then the distribution function of the fluid elements, f(t,x,u)du, 
can be determined independently of the solution of the kinetic equations 
for the energy and chemical species. In a compressible flow, determination 
of fdu from the existing solution of the momentum equations and solution 
of the kipetic equations become coupled. Since the existing solutions 
depend on the mean density variations only, an iterative scheme is expected 
to be satisfactory. 

The present report describes the initial phase of development of the 
approximate formalism. It consists of two portions. The first contains 
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the solution of the two nonequilibrium-degree kinetic equation ^ for the 
effective length scale and turbulence energy for a spatially homogeneous 
turbulence field with two characteristic length scales, wherein the source 
for one family of eddies exists. This solution is applied to the evalua- 
tion of the eddy diffusivity in the combustion chamber of an internal 
combustion engine. The result is compared with the other existing solu- 
tion. ^ This is carried out to demonstrate the feasibility of obtaining 
an effective length-scale equation within the context of the kinetic theory. 

The second portion of this report contains the formulation and partial 
solution of the compressible plane shear layer according to the approximate 
formalism discussed earlier. 


II. HOMOGENEOUS FIELD WITH TURBULENCE SOURCE 


Governing Equations 

We consider an incompressible, homogeneous turbulence field consisting 

of two characteristic length scales. At time zero, the turbulence energies 

associated with both length scales are given. For time greater than zero, 

we consider that a source for the eddies associated with one of the two 

scales exists in the field. The two coupled Langevin equations for the 

-*■ ■+ 

two degrees of freedom, U and V, given in References 5 and 6 are modified as, 



dt 


- B (U! - V?) + A. (t) 

c l i l 
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K 
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( 1 ) 
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- £ b V i 


B c (V i " U p + A i (t) 


+ 



( 2 ) 


where the U 1 -degree of freedom is considered to have the source for t > 0. 
Also, as it is evident in Eq. (1), the source strength, 2SUj, is considered 
to be proportional to U^. As we shall see subsequently, this particular 
form can imply, within the context of the present simplified analysis, the 
turbulence energy generation in proportion to the mean shear. 

All symbols are defined in the Nomenclature. The physical meaning of 
the Langevin equations has been discussed in the previous papers, and, 
therefore, it is not discussed here. 

The Fokker-Planck equation governing the joint distribution function, 
f 2 (t,x,U, , \^,)dU , dV , , can be constructed from Eqs. (1) and (2) in the same manner 
as that shown in Reference 5. The moment equations governing the three 
components of the turbulence energy then can be derived from the Fokker- 
Planck equation in the same manner as that shown in Reference 6 as, 


d<u kV 

dx 
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(5 + 


M)<u kV 


M<v kV 
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d<U kV 1/2 I 

-jr*- * ‘"kV 3<1 VV + (2 + “'W 


-[ 


7 + 3M - 


— Stj 3(t) ] 
‘W J 


<u kV 


( 5 ) 


where 


W ■ <u kV + 2< W + <v kV 


(6) 


In the above equations, all energy quantities are in dimensionless form 


as. 


<u kV 


<KK>/<KK> 

k k k k o 


<u kV 


<u i u*>/<w'w' > o 


( 7 ) 


<V kV ’ < ?k >/< Vi , o 


Also, 


M - A a /Ah 


S/P 


( 8 ) 


ao 


e ao ■ <W i H k^ (2A .> 


and 
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The three coupled equations, Eqs. (3), (4) and (5), can be readily 
integrated as soon as a set of initial conditions and the source function, 

A 

S(x), are specified. These values will be specified to simulate the 
physical problem found in the piston engines to be discussed in the sub- 
sequent subsection. 

From the solution of Eqs. (3), (4) and (5), the effective length 
scale, A e (t), can be constructed as (see References 5 and 6): 


A 

e 


<W. W. > 
k k 


[<u k u k > + <u u v v >] + m[<u 1 / v u > + <v w >] 


Vk 


k k 


Vk 


(9) 


In the present problem, A g represents the characteristic scale of the hypo- 
thetical family of the energy-containing eddies which would give the same 
time variation of the turbulence energy, <W k W k >, of the field as that 
caused by interaction of the two families of the eddies, U and V modes. 
(See Reference 5.) 


Eddy Diffusivity in Piston Engines 

As it was stated in the Introduction, one of the objectives of the 

present report is to demonstrate a feasible means of deducing a length-scale 

equation usable within the context of the kinetic theory approach. We 

shall employ the present two-nonequilibrium-degree equations, Eqs. (3) - 

(5), to analyze a portion of the problem found in a piston engine first 

(7) 


studied by Sirignano. 
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Consider the sketch of a piston engine cylinder shown In Figure 1. In 
order to analyze the combustion taking place in the cylinder, the turbulent 
mixing of the fuel and air must be described. This mixing in the cylinder 
is controlled by at least two major families of the characteristic eddies 
with considerably different characteristic scales. 

The fuel and air enter the cylinder turbulently through the valve. 
Therefore, their turbulence characteristics are of those dictated by the 
valve. Subsequently, however, larger scale eddies are generated by move- 
ment c: the piston. In the analysis of Reference 7, which is based on the 
mixing-length theory, the eddy dif fusivity, D'(t) , was approximated by the 
follow ng equation. 

D* (t) - A a I u^ j + A fa Uj exp[-Uj(t - t i ''/A b ] (10) 

where u* and u* denote the instantaneous piston velocity and the initial 
P I 

intake velocity of the gases through the valve, respectively. and A^ 
denote the scales of the eddies generated by the piston movement and the 
valve, respectively. 

The rationale for Eq. (10) is that the characteristics of the smaller 
eddies initially generated through the valve decay exponentially whereas 
the piston velocity itself represents the mean turbulence velocity of the 
larger eddies generated by the piston at any given time. 

In order to simulate this problem by the homogeneous solution of 
Eqs. (3), (4) and (5), we first consider that the V-mode of the turbulence 
represents the smaller scale eddies generated by the valve, and that the 
U-mode that of the larger scale eddies generated by the piston movement. 


Thus, 
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M - A /A. > 1 

A O 


(ID 


Then, ve assume that the piston movement engenders the eddies by creating 
a shear layer along the wall. (See Figure 1.) If we let A^ denote the 
scale which Is of the order of the shear layer, the mean velocity gradient 
will be of the order of (u^l/A^. It is well known ^ that the generation 
rate of turbulence by shear is proportional to the shear stress times the 
mean velocity gradient. Within the context of the present simulation this 


rate would be of the order of > 

k k 



For the present solution, we 


let 


S 



2 A 


A 


( 12 ) 


Then, 


S/0 


ao 


e "i7j 

< W. 


(13) 


and it can be readily shown from Eq. (3) that the rate of production of 

Note that there exists a rather basic difference between the assumptions 
of Eq. (10) made in Reference 7 and thone of the present equation. In Eq. 
(10), it is assumed that the piston velocity |u^| itself represents the 
instantaneous mean turbulence velocity generated by the piston movement. 

In the present assumption of Eqs. (3) - (5), and Eq. (13), the piston 
velocity, | | , is considered to be the velocity which is responsible for 
the mean shear which gives a rise to the turbulence source. The 
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instantaneous mean turbulence velocity generated by the piston movement is 

1 /2 

<U'U;> ' which will be found from the solution of the governing equations* 
Eqs. (3) - (6). 

From the solution of the present equations, variation of the effective 
length scale, A^, is obtained from Eq. (9). Then, again within the context 
of the present simulation, we write, 

D'(t) - VW ^> 1/2 (U) 

The eddy dif fuaivitiea given by Eqa. (10) and (14) are nondimensionalized as. 
Reference 7 

D(t) - I« p | + i < Vk >y 2 .,[>V kY f T ] (15) 

Present 

D(t) - l <W.W, > 1/2 (16) 

e k k 

where 

Up i ■ i^i/^y 2 , d(t) - D’(t)/A a < Wl ;w*>y 2 < i7 > 

In writing Eq. (15), we took cognizance of the fact that the initial 

1/2 

mean turbulence velocity of the intake gases, <V ]V^> 0 , is equal to the 

intake velocity u^. (See Equation 10.) 
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Solution of Eqs. (3), (4) and (S) 

The initial conditions of the present equations for <U^Uj t > and 
are chosen such that sons meaningful comparison can be made between the 
present solutions and Eq. (15) employed in Reference 7. It is pointed out 
at the outset, however, that a deatiled quantitative comparison between the 
two results is not possible. This is because there is no quantity In the 
expression for D'(t) of Reference 7 which corresponds to the present turbu- 
lence source strength, S. In Reference 7, as it has been already explained, 

1/2 

a direct assumption on the < U^U^> itself has been made as being equal 
to | | . In the present simulation, S was formulated by postulating a 
turbu) err. i generation mechanism (Figure 1), and specification of its value 
beyond that of the order of magnitude is meaningless. Nevertheless, as 
we shall see, the comparison will elucidate certain basic qualitative aspects 
of the multi-scale mixing phenomenon. 

For t - 0, we let, 


<V*V'> 1/2 - u! 

k k o I 


(18) 


and 


<u'u ' >1/2 

k k o 


<U’V’> 1/2 
k k o 


|u-| 


. 0 


in Eqs. (3) - (5). 

Figures 2 through 5 show the typical results. 


Discussion 

Figure 2 shows the contribution of the smaller eddies generated by the 
Intake valve, D v , to the diffusivity for the constant piston velocity, | ) • 
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S-B denotes the results of Reference 7 given by Eq. (15). Since it has been 
assumed in Reference 7 that the contribution of the smaller eddies decays 
exponentially, the broken lines merely show these exponential decays for 
the two given scale ratios, A^/A a , of 1/3 and 1/10, respectively. The 
present solutions, on the other hand, show that the contribution of the 
smaller eddies does not completely decay but it persists. This is due to 
the fact that the smaller eddies continuously receive energy from the 
larger eddies which are being generated by the piston movement. 

Figure 3 shows the diffusivity D generated by both the intake valve 
and piston. According to Eq. (10), D' monotonically decreases to the constant 
value, A a |u^|, as dissipates exponentially. In the present solution, 

D initially increases with t as the constant piston velocity begins to 
generate the eddies. However, D begins to decrease as D^-dissipation comes 
to exceed the piston generation of the larger eddies. For t of order one, 
the rate of generation of the larger-scale turbulence by u^ begins to be 
balanced by the turbulence dissipation of the smaller eddies, and D approaches 
an asymptotic value. As it has been alluded to, comparison, between Eq. (15) 
and the present solution, of the absolute levels of D(t) is rather meaningless. 

Figures 4 and 5 show D y (t) and D(t) for the sinusoidal variation of u^. 

Since the smaller eddies are assumed to be independent of u^ in Eq. (10), 

Figure 4 shows that D v decays exponentially in the S-B solution as it did 

for the constant u (see Figure 2). In the present solution, D somewhat 
p v 

follows the sinusoidal piston movement but with a time lag after the initial 
decay period. The main point here is that the contribution of the smaller 
eddies to the diffusivity persists for all t as the smaller eddies continuously 
receive energy from the larger eddies being generated by the piston movement. 
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Figure 5 shows the variations of D(t) for the sinusoidal variation of 
Up. In the S-B solution of Eq. (15), D y almost completely dissipates for 
t > n /^» ant * D becomes einusoldal as is u^. The present solution shows 
that there exists a substantial time lag between the diffusivity and the 
piston velocity since the u^Cx) determines the instantaneous turbulence 
source, S(x), but not the instantaneous turbulence energy. 


III. COMPRESSIBLE PLANE SHEAR LAYER 

As the initial development of the formalism discussed in Section I, 
which will at' > xse the existing solution of the mean momentum equation in 
the solution of the kinetic equations, the compressible plane shear layer 
is being analyzed. This section describes the analysis. The final numerical 
results are not as yet available. However, all basic features of the 
solution expected as well as the method and assumptions made are discussed 
herein. 


Existing Solution of the Momentum Equation 

There are several solutions of the mean flow properties available. For 
convenience of the present analysis, we shall employ the compressible solution 
obtainable via the Howarth transformation from the incompressible solution 
given in Reference 9. 

We first define the density-averaged ordinate, 


ty 


P- 


dy 




£ 

X 


o 


(19) 
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Then, we employ the Ting-Libby compressibility modification (Reference 10) 
of the Incompressible eddy dlffusivlty of momentum as, 


e 


0.014 bu 
00 



X 


( 20 ) 


With Eqs. (19) and (20), and with appropriate similarity assumptions, the 
governing mean momentum equation of the compressible plane shear layer 
based on the mixing-length theory becomes Identical to that for the incom- 
pressible flow whose solution is given in Reference 9. The momentum equa- 
tion transformed on the similarity plane is then, 




0 


( 21 ) 


with the boundary conditions, 


dF 

dC 


(-) 


1 + 


u 

00 
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—GO 


+ U 


dF 

d£ 


(_oo) 


- 1 


u„ - U _ 

go _ao 


u + u 
00 *00 


( 22 ) 


From the solution of Eq. (21), the mean velocity and shear stress profiles 
can be established as. 


u + u 
00 *00 


<U> 
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l + 


u - u 


u + u 
00 *00 


erf i 
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-i £ ( 


(23) 


<v> * 


(24) 
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t 


<UV> - - (0.014) bo 2 f (25) 

\ 00 / d£ 

Reference 9 suggests that the values of the constants o and b be 13.5 
and 0.098, respectively. The functions, F and dF/d£» are obtained from 
Reference 9. 


Kinetic Equation for Compressible Flow 

There are substantial experimental indications (see References 11, 12 
and 13), which seem to show that combustion does not alter the fundamental 
structure of the turbulence field. In the absence of definite contrary 
experimental data, we shall assume that the high-turbulence-Reynolds-number 
energy cascading structure of the turbulence field upon which the incom- 
pressible kinetic theory has been built remains valid. Then, the Vinetic 
equation can be constructed to include the density variation in basically 
the same manner as the incompressible kinetic equation of Reference 1. 

This compressible equation is. 


(< v 
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3pzf 

3x k 


3<u > 
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3pzf 

3U 
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3U, 


( f zU k ) 


+ 


<U U > -2- 1 

m m 3 fz 

3 3U, 3U. 

k k J 




tuT (r k £ «> 


- fpw • 0 


(26) 


The structure of the turbulence field enters into the kinetic equation 
through the eddy-interaction terms contained in the square brackets, and it 
is seen that these terms are identical to those derived for the Incompressible 
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kinetic equation since it has been assumed that the fundamental structure 
of the turbulence field is unaltered by the compressibility. 

r k , in Eq. (26), represents the direct turbulence generation by the 
compressibility. This term should engender, among other terms, the corre- 


lation 


< P 


t— - > in the second order moment of Eq. 


(26). 


It is our 


present feeling that this term is not important unless there exist shocks 
in the flow field (Reference 14). Therefore, for the present initial 
analysis of flows without shock?, we shall let 


will be analyzed in detail in a future paper. 

Once the solution of Eq. (26) is obtained for a given flow problem, all 
one-point moments of engineering interest can be readily constructed as, 
for instance, 

<u i> ■ j fu^dU 

(27) 

<pU 1 U j > - | pfct.U)!^ f (x,U)dU 

Before Eq. (26) can be solved, a specific relationship must be found 
between the density, temperature, and the pressure. 

There are several plausible relationships that can be employed to relate 
the fluctuating and mean components of p, t, and P. For the present study, 
we shall try the relationship suggested in Reference 15 for ideal gases, 
which is, 
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t - <t> 

<t> 


P ~ <P> 

<p> 


(28) 


where t - <t> and p - <p> denote the fluctuating components of the temperature 
and density, respectively. Equation (28) can be rewritten as, 


<p> 



(29) 


Solution of Compressible Plane Shear Layer 

According to the approximate concept set forth in Section I, we shall 
make use of the existing solutions given by Eqs. (23), (24), and (25). The 
particular problem to be analyzed is the plane shear layer shown in Figure 6. 
In addition to the mean velocity difference, u and u , it is considered 
that the two streams are at different mean temperatures of t w and t^, respec- 
tively. We seek the solution of the kinetic equation, Eq. (26), for the 
temperature and density profiles, as well as the other turbulence quantities, 
with the aid of Eqs. (23) - (25). 

As an initial try, we assume that the probability distribution function 
of the fluid elements, f, and those of the temperature and density, tf and 
pf, can be approximated by bimodal functions. Thus, we write, 

f - f l + f 2 (30) 


where 
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for V > 0 


for V < 0 


(31) 


(32) 


with f^ ■ 0 for V ^ 0 and fj ■ 0 for V ^ 0, as it has been done earlier. 

(See References 1-5.) 

Using the assumed bimodal form of the distribution function of Eq. (30), 
we derive, 


<u> 



+ U o2> 


(33) 


< v> - _L_ (/T“ - /el ) 
^7 1 2 


(34) 


<BV> ' 7 ^ <u » 2 ' u ° 2> (/i 7 + /I 7 ) <35) 

- <u kV ■ + + (i-k) < E i + V 

- + h m> 


The first step of the present approximate solution is to determine the 
functions u q ^( 0» 0^(0* E^(Oi and Ej^) by matching Eqs. (33) - (36) with 
the corresponding moments obtained from the existing solution of the momentum 
equation. Equations (23), (7.4) and (25) give <u>, <v>, and <UV>, as functions 
of £. However, the first order mixing-length theory, from which these 
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functions have been obtained, does not give the turbulence energy, 
We shall assume that 


<U kV " ~ 5 <UV> (37) 

Equations (23), (24), (25) and (37), show that these four mean quantities 
are given provided that the density ratio P 00 / < P > is known. This is to be 
expected since in a compressible flow the momentum and energy equations are 
coupled. 

The present proposed method of solution is as follows. For a given problem 
with u /u and t /t ("p /p ), the density ratio, <p>/p , is assumed. Equa- 

—OP OP —OP 00 OQ —OP 00 

tions (23) - (25) and (37) then give the four moments, <u>, <v>, <UV>, and 
< U k U k > as functions of £. With these moments given, Eqs. (33) - (36) consti- 
tute a set of four nonlinear algebraic equations for the four unknown functions, 
u o l(0, u 0 2^* *^(0» and ^ 2^0 • Solution of these equations for the four 
functions results in the distribution function f through Eq. (30). With the 
f known, the kinetic equation, Eq. (26), is solved for the distribution 
function of the density, and, therefore, of the temperature through Eq. (29). 

The mean density profile is then evaluated as, 




1_ 

P« 


J 


p(t,V) f(t,v)dv 


(38) 


With the use cf the < P>/p oo evaluated by the above equation, the moments, 
<u>, <v>, <UV>, and are recalculated from Eqs. (23) - (25) and (37). 

Solution of Eq. (26) is then repeated. This iteration continues until a 
satisfactory convergence on the density is attained. 

We shall now discuss the solution of the kinetic equation, Eq. (26). 
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Commensurate with the blmodal assumption of f (see Gq. 30), Eq. (26) Is 
solved also by an approximate blmodal method. As It has been done earlier 
(References 1 - 5) , we let, 


t(C,$) - t 1 (c,v) + t 2 U,v) 


(39) 


P<£,U) - Pl tt,v) + P 2 (*t,V) 


(40) 


where t^ and p^ are zero for V <_ 0, and t 2 and p 2 are zero for V >_ 0. 

Now, we generate two moment equations from Eq. (26) in the following 
manner. We first let z ■ t and integrate Eq, (26) termwise with respect 
to the velocity space. This produces the thermal energy conservation equa- 
tion. Next, we multiply Eq. (26) through by V, with z * t, f id Integrate 
the resulting equation termwise with respect to the velocity space. This 
results in the turbulent thermal-energy transport flux equation. In both 
of the above equations, we set - u - 0. The resulting moment equations 
are as follows where the appropriate similarity transformation commensurate 
with the momentum equation, Eq. (21), has been performed. 
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where 
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OV> - —I: (e^O,) (^+^> (48) 

2/67 1 L 1 L 


9 1 " V*- . 0 2 » t 2 /t ai (49) 


In Eqs, (41) through (48), the functions, <u>, u u q2 , <v>, 
and , represent the dimensionless quantities which have been normalized 
by u n . Also, p represents a dimensionless quantity which has been 
normalized by p m . The same symbols as for the previous dimensional quanti- 
ties have been employed for convenience. 
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Solutions 

A typical solution of Eqs. (33) - (36) for the functions, u^, u^* 
t^E^and /EJ comprising the distribution function of the fluid elements f 
Is shown In Figure 7, for a given density ratio <p>/p a) . 

As it was mentioned earlier, solution of the kinetic equations, Eqs. 
(41) and (42), has not yet been completed. Different numerical techniques 
to accelerate the convergence, with respect to <p>/p oo , are being tried. 
Detailed numerical analysis of these equations will constitute a future 
paper. 


IV. CONCLUDING REMARKS 

As an initial development of a formalism which would substantially 
simplify application of the kinetic theory concept to the turbulent chemically 
reacting flow problems of engineering interest, the following two problems 
have been analysed. 

For the purpose of eventually deducing a length-scale equation 
commensurate with the kinetic theory, the two-nonequilibriura degree kinetic 
equation developed earlier has been solved for a homogeneous turbulence 
field with two characteristic families of the energy-containing eddies 
wherein a turbulence source exists for one of the two families. Solutions 
have been obtained which would simulate the turbulence which exists in 
the combustion chamber of an internal combustion engine. There, the family 
of the smaller eddies generated through the intake valve decays while the 
larger-scale eddies are continuously generated by the subsequent piston 


movement . 
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A comparison was made batween ths eddy diffusivltics deduced from the 
present solution and those obtained elsewhere by the mixing-length theory. 
Among other things, the present solution shows that the contribution of 
the smaller eddies to the dlffusivlty does not quickly dissipate as was 
predicted by the existing mixing-length theory. Through the interaction 
of the eddies, It was found that the smaller eddies continuously receive 
energy from the larger eddies In such a way that the effect of the smaller 
eddies on the dlffusivlty persists. 

The governing equations deduced relate the effective length scale, A g , 
to the turbulence energy, turbulence-energy source strength, and the two 
given length scales of the two families of the eddies. It seems very 
feasible that such equation could be solved in conjunction with the 
simplified formalism being developed. 

The fundamental concepts of the simplified formalism being proposed have 
been discussed in the text. These concepts are being applied to the solution 
of a compressible plane shear layer. A simple, bimodal approximation has 
been employed to deduce the distribution function from the existing, mixing- 
length solution of the momentum equation, A governing set of the moment 
equations for the temperature and density fields has been constructed 
according to the kinetic theory by the use of the distribution function. 

This set of the equations is coupled with the existing solution of the 
momentum equation through the density. 

Numerical solution has not yet been completed. Different approximations 
will be tested in the future for the deduction of the distribution function 
of the fluid-elements from the existing solution of the momentum equation. 
This seems to be one of the key points which may dictate the outcome of the 
proposed formalism. 
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NOMENCLATURE 


A t (t) 


b 

U(t> 

D’(t) 

D v (t) 

E l» E 2’ 

F(0 

f 


f l» f 2 


f 2 (t, x , U\ V') 


..l* S.: 


Random acceleration of fluid element by the small equilibrium 
ed< lee. 

constant for eddy dlffuslvity 

dimensionless dlffuslvity defined by Eq. (17). 

eddy dlffuslvity. 

contribution to D(t) by the valve-generated smaller eddies, 
functions defining f^ and f^ respectively, 
dimensionless stream function for plane shear layer, 
distribution function of fluid elements, 
functions defined by Eq. (30). 

one-point joint distribution function for U'and V' modes, 
forces due to molecular viscosity. 


M 

P 

P 

S 

A 

s 

t 

C I 

u 

u\ u\ 

U.U 4 


A /A, 
a b 

flucluatlng portion of pressure. 

Instantaneous pressure 

turbulence source strength due to piston movement. 

S/6 .o 

time, or temperature (Sec. Ill) 
time at which value opens 

x-component of the velocity relative to the mean velocity in 
Sec. in. 

U’- mode velocity. Contribution of the family of larger eddies 
‘ to W\ 


‘W* <u kY J 

<v kY’ <M kV > 


<u.u.> 

k k 


relative velocity in sec. III. 

dimensionless energies 
defined by Eqs. (7). 
Turbulence energy in sec. III. 





i 


U 


u, u A 


u' 


U ol‘ U o2 

V 

v' V* 
v, v A 


v 

W 

w\ W| 

-*■ 

x, x 1 

X, y, z 
z 
6 


V *b 


ft , A. 
a b 


e 

Ae 
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x-component of the absolute velocity. 

instantaneous absolute velocity in sec. III. 

mean turbulence velocity of the value-generated eddies. 

piston velocity 
1/2 

u' <WW> ’ 
p k k o 

functions defined by Eqs. (31) and (32). 

y-component of the velocity relative to the mean velocity. 

V'-igode velocity. Contribution of the family of smaller eddies 
to V. 

relative velocity in sec. III. 
y-component of the absolute velocity. 

z -component of the velocity relative to the mean velocity in sec. III. 

velocity relative to the mean velocity. 

postion vector. 

cartesion components of x 

scalar quantity. 

<U k U k >JS/2A ln S6C * 111 * 

<W ’k W V^ 2A a and <W ’k W ’k >55/2A b re8 P ectivel y* 


dimensionless quantity defined by Eq. (43) 
direct effect of compressibility on turbulence 
eddy diffusivity 
t/t w 

characteristic length scale of plane shear layer (*x/c). 

characteristic scales of U 1 and V’ modes respectively. 

effective scale 

A /A 
e a 

similarity variable defined by Eq. (19) 
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p density 

o constant 

t 38 t/2 

ao 

to source for z due to chemical reaction 

< > ensemble average 

Subscripts 

o initial 

«, -» averaged quantities at y ■ ■ and y « 


respectively. 



tmnM zmrm jjovj 


VNianmi] 
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v 


e/t- B v/ q V J0 * (1) Aq 









(T) for A b /A a =1/10 



(T) for A b /A a =1/3 
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Fig 5. 


Diffusivity for sinusoidal piston velocity. 








